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Overview  of  the  research  objectives 

The  objective  of  this  research  effort  was  to  develop  the  expertise  and  the  experimental 
research  capabilities  for  precision  laser  Doppler  measurements  of  surface  vibrations 
relevant  in  stractural  acoustics.  Over  the  three  years  covered  by  the  AASERT  grant  (07/94 
to  08/97),  the  main  research  effort  concentrated  on  building  a  laser  system  that  could  be 
used  to  determine  the  two  complex  elastic  moduli  (bulk  and  shear)  of  a  voided  polymer 
within  known  and  acceptable  error  bounds. 


Relevance  to  the  Navy 

Effective  ship  hull  treatment  depends  in  part  on  a  good  understanding  of  the  acoustic 
properties  of  the  viscoelastic  layer  attached  to  the  hull.  Measuring  the  shear  and  bulk 
dynamic  moduli  of  a  viscoelastic  sample  is  therefore  critical  to  predict  correctly  acoustic 
performance.  Two  of  the  most  common  experimental  techniques  used  to  measure  dynanuc 
moduli  are  the  resonant  technique  (Madigosky  and  Lee)  and  the  DMTA  (Dynamic 
mechanical  testing  apparatus),  both  of  which  are  very  useful  but  have  some  intrinsic 
limitations.  The  proposed  laser-based  experimental/numerical  technique  is  an  alternative, 
noncontact,  method  capable  of  measuring  simultaneously  both  shear  and  bulk  moduli.,  an 
important  advantage  since  manufacturing  variabilities  and  hysteretic  material  response  are 
common  sources  of  errors. 


Technical  Approach 

In  the  first  phase  of  the  research  program,  a  set  of  five  independent  laser  Doppler 
interferometers  was  built,  calibrated,  and  tested  for  measurements  of  stmctural  vibrations  in 
the  0.3  -  5  kHz  range.  The  design  relies  on  optical  fibers  so  that  the  laser  probe  heads  can 
be  placed  relatively  far  from  the  optical  bench.  The  optical  probe  heads  were  designed  to  be 
very  compact  and  capable  of  measuring  either  out-of-plane  or  in-plane  surface  motion. 

(Two  photographs  of  the  experimental  setups  are  available  on  the  web  site 
http://www.me.gatech.edu/yves.berthelot/bothphotos.html) 

The  second  phase  of  the  research  was  to  use  the  system  to  measure  the  frequency- 
dependent  bulk  and  shear  moduli  of  some  voided  polymers.  A  sample  of  arbitrary  shape 
(typically,  a  rectangular  block  of  r’x2”x3”)  is  excited  harmonically  at  its  base  (by  a  shaker 
or  a  piezoelectric  stack).  The  surface  dynamics  of  the  sample  is  non-invasively  measured  at 
several  points  (in-plane  and  out-of  plane  surface  velocities  -  amplitude  and  phase)  by  a  set 
of  four  independent  fiberoptic  laser  Doppler  interferometers  and  referenced  (amplitude  and 
phase)  to  its  base  motion.  The  material  parameters  are  determined  from  a  le^t-square  fit 
between  the  experimental  data  and  numerical  predictions  obtained  from  a  finite  element 
code  in  which  the  frequency-dependent  elastic  moduli  are  the  adjustable  parameters. 


Major  Results 

1)  Hardware 

2.1:  Five  probe  laser  system 

A  five  probe  fiberoptic  interferometer  system  has  been  designed  and  calibrated  for 
laboratory  measurements  of  surface  vibrations  at  ambient  conditions  has  been  built  and 
calibrated.  (See  Willis  et  al.  JASA,  102  (5),  Nov  1997).  See  also  Figure  1.  Typical 
performances  and  characteristics  are  shown  in  Table  1. 


In-plane 

Out-of-plane 

Noise  floor  (w/o  ave.) 
Noise  floor  (w/64  ave) 
Maximum  velocity 
Maximum  displacement 
Sensitivity 

Bandwidth 

65  ^m/s  -  rms 

5  |im/s  -  rms 

35mm/s  -  peak 
18p.m-pe^ 

23.7  (|xm/s)/mv 

DC  - 10  kHz 

35  |j.m/s  -  rms 

3.5  pm/s  -  rms 

17  mm/s  -  peak 

125  pm -peak 

11.7  (pm/s)/mv 

DC  - 10  kHz 

Table  1:  Typical  pe/formance  and  characteristics  of  the  interferometers. 


1.2:  Two  probe  scanning  system 

A  two  probe  scanning  system  is  being  completed  for  measurements  inside  a  small 
pressure/temperature  chamber.  (0-500  psi ,  6-40  C).Direct  measurements  in  a  pressure  - 
temperature  chamber  will  alleviate  some  of  the  problems  currently  encountered  with  the  use 
of  time-temperature  superposition  principle  and  associated  shifts. 

1.3:  Single  probe  system 

A  separate  single  proble  system  (either  in-plane  or  out-of-plane)  for  independent 
calibrations  and  testing.  This  system  can  be  also  used  in  the  Mhz  range  for  ultrasonics. 


2)  Software 

2-1:  Automated  data  acquisition  (Labview)  and  Processing  (Matlab) 

The  5  laser  probe  system  has  been  automated  for  data  acquisition  and  processing.  The 
system  is  capable  of  acquiring  7  signals  (from  two  digital  oscilloscopes),  including  2 
reference  signals,  64  times  at  a  given  frequency.  Each  signal  is  analyzed  with  a  DFT  to 
measure  the  complex  amplitude  at  the  frequency  of  interest  relative  to  the  complex 
amplitude  of  the  base  excitation. 

2-2:  Finite  element  code  (Fortran) 

A  finite  element  code  is  used  to  model  the  dynamic  behavior  of  the  sample.  The  code  uses 
825  nodes,  and  36  quadratic  elements.  The  boundary  conditions  are  such  that  the  all  faces 


of  the  sample  are  free  except  for  the  bottom  face  which  is  driven  harmoriically  and 
uniformly.  The  model  assumes  a  homogeneous  and  uniform  material  with  Young  modulus 

E,  loss  factor  ti,  and  Poisson  ratio,  v  . 

2.3:  Inversion  Code  (Fortran) 

The  search  for  the  3  parameters  (E,  T|,  v)  that  best  describe  the  measured  data  is  obtained 
by  a  three-dimensional  direction-set  method  (Powell’s  method),  a  robust,  classical  method 
in  optimization  theory.  The  function  that  is  minimized  is  the  mean-square  difference 
between  the  measured  values  (amplitudes  and  phases)  of  surface  displaceinents  relative  to 
the  base  and  the  values  predicted  by  the  FEM  code.  (See  block  diagram  Figure  2). 

For  more  details,  see  Willis  et  al.,  JASA  102  (5),  Nov  1997. 


3'>  Experimental  results 

3-1:  Frequency  dependent  bulk  and  shear  moduli  and  complex  sound  speeds 

The  frequency-dependent,  complex  sound  speeds  (shear  and  longitudinal)  were 
simultaneously  measured  in  the  same  sample  in  the  0.5  -  2.5  kHz  range.  (The  frequency 
range  could  be  lowered  with  a  larger  sample).  The  results  are  shown  in  Figure  3.  (solid 
line).  Note  that  the  vertical  scale  is  greatly  expanded  around  the  mean  value,  thus  indicating 
the  great  precision  of  the  technique. 

3-2:  Error  analysis 


•  Intrinsic  variability  (repeatability):  A  statistical  analysis  of  the  data,  over  64  sets  of 
data,  was  performed  to  measure  the  variance  of  the  data  obtained  with  the  same  sample 
and  the  same  measurement  apparatus.  The  resulting  variance  in  the  complex  sound 
speeds  was  estimated,  see  Figure  1  (dashed  lines).  The  dashed  lines  correspond  to  the 
sound  speeds  evaluated  with  the  mean  surface  velocities  ±  1  standard  deviation.  This 
shows  that  the  error  in  amplitudes  and  phases  of  the  measured  surface  velocities 
translate  into  variation  of  sound  speeds  less  than  1%. 


Relative  error  (*)(l  .0-2.5  kHz) 

real  part 

imaginary 

Longitudinal  sound  speed 

Shear  sound  speed 

<1%  (1%) 
<1%  (1%) 

<1%  (3%) 

<1%  (3%) 

(*)  values  indicate  the  typical  average  error  over  the  frequency  range.  Values  in 
parentheses  indicate  the  maximum  error  in  that  frequency  range. 


•  True  variability  between  measurements  (same  sample):  Tested  the  effect  of  the 
variability  of  bonding  the  sample  to  the  base,  i.e.,  repeated  the  measurements  with  the 
same  sample,  before  debonding  and  after  rebonding  to  the  same  base.  In  each  case,  the 
base  was  excited  with  a  different  drive:  a  Ling  shaker  ^d  a  piezoelectric  stack.  The 
relative  errors  measured  with  the  same  sample  tested  with  the  two  drives  (and  with 
debonding  and  rebonding  to  the  base)  are  shown  in  the  table  below. 


Relative  error  (*)(i. 0-2.5  kHz) 

real  part 

imaginary 

Longitudinal  sound  speed 

Shear  sound  speed 

3%  (4%) 

3%  (4%) 

10%  (13%) 

12%  (20%) 

(*)  values  indicate  the  typical  average  error  over  the  frequency  range.  Values  in 
parentheses  indicate  the  maximum  error  in  that  frequency  range. 


•  True  variability  between  measurements  (different  samples):  Tested  the  method 
with  two  samples  taken  from  the  same  piece  of  material,  both  with  the  same 
excitation  (Ling  shaker). 


Relative  error  (*)(l. 0-2.5  kHz) 

real  part 

imaginary 

Longitudinal  sound  speed 

Shear  sound  speed 

4%  (5%) 

4%  (5%) 

10%  (16%) 
10%  (16%) 

(*)  values  indicate  the  typical  average  error  over  the  frequency  range.  Values  in 
parentheses  indicate  the  maximum  error  in  that  frequency  range. 

Comaprison  of  the  results  with  independent  (in-situ)  measurements  made  by  Jarzynski 
et  al.  on  a  large  chunk  of  the  same  material.  Also,  compared  the  results  with  Kemer’s 
model  in  which  the  matrix  material  properties  were  calculated  using  experimental  data 
obtained  with  the  DMTA  apparatus  (results  provided  by  Walt  Madigosky). 

results:  (a)  In  situ  results:  Sound  speeds  are  in  reasonable  agreement  with  our 
measurements,  but  at  present,  the  in-situ  results  are  only  accurate  within  10  to 
20%  and  only  the  real  partof  sound  speeds  has  been  estimated. 

(b)  Kemer’s  model.  The  comparison  with  Kemer’s  simple  model  of  voided 
polymers  is  not  really  meaningful  because  we  don’t  know  the  material  properties 
of  the  host  material.  (We  had  to  estimate  the  properties  from  measurements  made 
with  the  DMTA,  measurements  provided  by  W.  Madigosky).  Nevertheless,  the 
comparison  between  measured  values  of  sound  speeds  and  predicted  values  was 
remarkably  close. 

Spatial  uniformity  of  the  base  motion  (0.5  to  3  kHz).  We  compared  the  uniformity  of 
several  bases  (inlcuding  a  wedge),  and  excitation  by  the  Ling  shaker,  the  B&K  4810 
shaker,  and  several  piezoelectric  stacks. 

result:  The  piezoelectric  stack  (Sensor  Technology)  performs  best.  A  simple 
aluminum  base  to  support  the  sample  works  best. 

Base  assembly:  We  tested  for  the  best  combination  of  bondings  between  the 
piezoelectric  stack,  the  back  mass,  the  base,  and  the  sample. 

result:  best  results  were  obtained  with  the  piezoelectric  stack  bonded  to  a  massive 
lead  brick  using  uncured  butyl  rubber;  the  aluminum  base  supporting  the  sample 
was  attached  to  the  stack  with  rubber  cement;  and  the  sampled  was  super-glued  to 
the  base.  With  this  arrangement,  the  nonuniformity  of  the  base  motion  was  less 
than  2%  below  2  kHz  and  less  than  4%  at  3  kHz.  (Results  with  the  Ling  shaker 
were  worse.) 


•  Reflective  tape.  We  assessed  the  effects  of  using  small  reflective  tapes  or  small  spots  of 
white  paint  on  the  sample,  for  either  in-plane  or  out-of-plane  measurements. 

results:  the  effect  appears  negligible. 

•  Sensitivity  analysis  on  the  position  of  each  laser  probe'.  We  performed  numerical 
experiments  (with  measured  data)  to  assess  the  relative  contributions  of  the 
measurements  made  by  each  laser  probe  on  the  resulting  material  properties. 

•  Sensitivity  analysis  of  the  numerical  code'.  We  performed  numerical  experiments  to  test 
the  dependence  of  the  results  on  (a)  the  initial  guesses  of  material  properties  (Young’s 
modulus,  loss  factor,  and  Poisson  ratio);  and  (b)  the  mesh  size  in  the  finite  element 
code. 

results:  the  material  properties  are  independent  of  the  initial  guesses.  They  remain  also 
identical  if  the  mesh  size  is  increased,  thus  indicating  that  the  method  is  very  robust  and 
does  not  converge  towards  false  minima. 
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Figure  1:  Experimental  setup 


@  each  frequency, 
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Figure  2:  Block  diagram  of  the  experimental  /  numerical  method 


Frequency-dependence  of  the  complex  sound  speeds 


solid  line:  sound  speeds  evaluated  with  the  mean  value  of  the  surface  velocities 

(64  acquisitions  at  each  of  the  5  surface  points  measured) 

dotted  lines:  sound  speeds  evaluated  with  the  mean  value  ±  1  standard  deviation 

of  the  sMace  velocities  (64  acquisitions  at  each  of  the  5  surface  points) 
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Figure  S 


Variability  between  2  samples  from  the  same  batch. 
Frequency-dependence  of  the  complex  sound  speeds 


solid  line;  sample  #1  -  dashed  line;  sample  #2 
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Figure.  4 
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